A light-front coupled cluster metho 



John R. Hiller and Sophia S. Chabysheva 

Department of Physics 
University of Minnesota- Duluth 
Duluth, Minnesota 55812 
(Dated: October 26, 2011) 

Abstract 

A new method for the nonperturbative solution of quantum field theories is described. The 
method adapts the exponential-operator technique of the standard many-body coupled-cluster 
method to the Fock-space eigenvalue problem for light-front Hamiltonians. This leads to an effective 
eigenvalue problem in the valence Fock sector and a set of nonlinear integral equations for the 
functions that define the exponential operator. The approach avoids at least some of the difficulties 
associated with the Fock-space truncation usually used. 
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I. INTRODUCTION 



Nonperturbative Hamiltonian methods that rely on Fock-space truncation are plagued by 
uncanceled divergences. The truncation removes contributions from couplings to higher Fock 
sectors that would otherwise cancel against contributions that are kept. When these con- 
tributions separately diverge with respect to an infinite-regulator-scale limit, an uncanceled 
divergence is born. For example, the Ward identity in QED requires contributions that differ 
in the number of photons required, as illustrated in Fig. [H a truncation in photon number 
will break the Ward identity [l| and introduce an uncanceled divergence. 

One approach to the handling of uncanceled divergences is to not take the infinite- 
regulator limit and instead seek a range of regulator scales over which physical quantities 
are slowly varying [2|. In this range, one has a compromise between errors due to trun- 
cation, which increase as the regulator scale is increased, and errors due to the presence 
of the regulator, which decrease as the scale is increased. This approach has been used in 
nonperturbative calculations of the anomalous magnetic moment of the electron regulated 
by inclusion of Pauli-Villars (PV) fields gQ. 

An alternative approach is sector-dependent renormalization where the bare param- 
eters are allowed to depend on the Fock sector. This absorbs the uncanceled divergences 
into redefinitions of the coupling constants. Calculations of the anomalous moment have 
been done this way jH, 6|. However, the Fock-state wave functions become ill-defined without 
proper normalization 7f , and the bare coupling is driven to zero at finite regulator scales jsf. 
In fact, the sector-dependent charge renormalization is not ordinary charge renormalization; 
it occurs in truncated QED without the presence of fermion loops and is the result of wave 
function renormalization in the broken Ward identity. 

To avoid these difficulties we consider a new method, the light-front coupled-cluster 
(LFCC) method @], that does not involve Fock-space truncation. It is instead based on the 
exponential-operator technique of the standard many-body coupled-cluster (CC) method ^ 
frequently employed in nuclear and chemical physics |lO|. We construct an eigenstate of 
a light- front Hamiltonian as \/Ze^\(f)), where |0) is a valence state with a minimum num- 
ber of constituents, T is an operator that increases particle number without violating any 
conservation laws, and is a normalization factor. The truncation that is made is in T 
but not in the exponentiation, so that all relevant Fock states are retained. This infinite 
expansion in Fock states necessitates some care in the computation of matrix elements, to 
avoid doing an infinite sum, but the techniques of the standard CC method can be used, 
with some extensions, to do the computations. 

In the standard CC method 10) is a product of single-particle states and terms in T 




FIG. 1. Graphs that contribute to the Ward identity at lowest order. The dashed vertical lines 
mark intermediate states. Only the first graph survives a truncation to Fock states with only one 
photon. 
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annihilate states in |0) and create excited states, to build in correlations. The operator T 
is then truncated at some number of excitations, usually no more than four. The focus is 
on finding the solution in the valence sector with a large but fixed number of particles. In 
contrast, in the LFCC method, the focus is on finding a solution in the infinite Fock space 
that has the same quantum numbers as the valence state The number of particles in 
the valence state not large, and how this state is to be determined is left unspecified. 



Some applications to field theory of the CC method have been previously considered [11 
for Fock-state expansions in equal-time quantization. The focus is on the non-trivial vacuum 
structure, which is avoided in light-front quantization, and particle states are then built on 
the vacuum. There was some success in the analysis of particularly of symmetry- 

breaking effects. 

For the LFCC method, light-front quantization is crucial, not only for the simple vacuum 
but also because the internal motion of a coupled system can be factored from the external 
motion. Of course, this separation can be done in the nonrelativistic CC method but not 
for a relativistic method in equal-time quantization. Obviously, we then use light-cone 



coordinates [12|, Il3[, which we define as = t + z for time, x_ = {x~,x±) for space, with 
x~ = t — z and x± = {x,y). The light-cone energy is p~ = E — pz, and the momentum is 
P = ip^iVi^i with p+ = E + Pz and pj_ = {px,Py)- These bring the mass-shell condition 

= to the form p~ = . The Hamiltonian eigenvalue problem is 

/Lf 2 _|_ p2 

v-m= ^ m- (1-1) 

In the next Section, we give a more detailed description of the LFCC method. This 



is followed in Sec. IIIII by a formulation of light-front QED in an arbitrary gauge 15|. In 
Sec. IIVI the two formulations are applied to a calculation of the electron's anomalous mo- 
ment. A brief summary is given in Sec. |Vl Additional details can be found elsewhere in 



these proceedings [14 . 



II. LFCC METHOD 

We wish to solve the fundamental eigenvalue problem (11. ip and compute physical quan- 
tities from matrix elements between eigenstates. The starting point is to write an eigenstate 
as = \fZe^\(j)). The factor Z maintains the chosen normalization, = S{P^ ~ E)- 

The valence state |0) is also normalized: = 5{P^ — P_). The operator T contains 

terms that only increase particle number, such that it conserves Jz, light-front momentum 
-P, charge, and any other conserved quantum number. Because is always positive, T must 
include annihilation, and, therefore, powers of T include contractions. These contractions 
will introduce loop corrections. 

To solve the eigenvalue problem, we construct an effective Hamiltonian V— = e~^V~e^ 
and let project onto the valence Fock sector. With these definitions, we arrive at the 
coupled system 

/Lf 2 _|_ p2 

P.V-\(j,)= p_, ^ 10), (1 - P.)Pi0) = 0. (2.1) 

The second equation determines the operator T, which then fully defines the effective Hamil- 
tonian in the first equation. The two equations are essentially a coupled set of integral 
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equations for the valence state and the functions that define the terms in T. In general, 
they must be solved by numerical methods. 

The expectation value for an operator O is given by 



Direct computation requires infinite sums. Following the techniques of the CC method [10 
we instead define O = e~^Oe^ and 



e e 



le^^e^l 



(2.3) 



We then find that (O) = (^|O|0) and 



im = (<l^'\7J^^^\<f>) = ^(P^-P)■ (2-4) 
The effective operator O can be computed from its Baker-Hausdorff expansion 

= d + [d,T] + ^[[d,T],T] + --- . (2.5) 
The bra {ip\ is a left eigenvector of as can be seen from the following: 

We still have an infinite system of equations. To arrive at a finite set, we truncate T at a 
fixed increase in particle count and truncate the projection (1 — Py) to a set of Fock sectors 
above the valence sector that yield enough equations to solve for the functions in T. This 
leaves a finite system of nonlinear equations for the functions in T and in the valence state. 
The effective Hamiltonian V~ = e~'^V~e^ has only a finite number of terms in its Baker- 
Hausdorff expansion, = V~ + T] + ■ • ■ and is computed from these commutators of 
T and V~ . Similarly, the expansion of an effective operator O has only a finite number of 
terms. The left-hand eigenstate \ip) is also truncated to a consistent set of Fock sectors, no 
more than found in T\(h). 



III. QED IN AN ARBITRARY COVARIANT GAUGE 

In order to facilitate a check of gauge invariance for physical observables, we quantize 
light-front QED in an arbitrary covariant gauge 15]. The PV-regulated Lagrangian is 



1=0 



1=0 

(3.1) 



A" 
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where the fermion and photon fields V'o and are the physical fields and indices of 1 and 
2 indicate PV fields. The interaction term is built from null combinations 



(3.2) 



t=0 i=0 

The coupling coefficients satisfy the constraints 

2 2 



eo = i, j2^-m = o, /3o = i, 5^(-im = o, 



(3.3) 



i=0 



i=0 



with ^2 and P2 fixed by chiral symmetry restoration 15, 1^ and zero photon mass 17 . 

To quantize, we apply a light-front analog of Stueckelberg quantization [18]. Consider 
the Lagrangian of a free massive vector field: £ = —^F"^ + ^fiA^ — \C{d ■ A)"^. The field 
equation is (□ + /i^)A/^ — (1 — Qd^{d ■ A) = 0, and the light-front Hamiltonian density is 



n = H\c=i + -{l-0{d- A){d- A- 2d-A+ - 2d± ■ Aj 



(3.4) 



with n\c=i = lT.l=o^^[{^±A^'y + ^'^{A'^y] and = (-1,1,1,1). The field equation i 
satisfied by 



Af,{x) 



dk 



. A=l 



^eW(A;) aA(fc)e-'*^-" + ai(fc)e''=" +ef(fc) ao(fc)e-''=-" + aS(A;) 



(3-5) 

with = /i for A = 1,2,3, but /iq = /i = [^l\fC,-, k = and k~ = (/c^ + iJ?)/k^. The 
polarization vectors 

e(i'2)(A;) = (0,2ei,2-^±A+,ei,2), e^'\k) = -{{kl - fx^) /k^ , k^ ,k^), e^''\k) = k/ fx, (3.6) 

/i 

satisfy k ■ e^'^-' = and e'-'^-' ■ e^'^'-' = — 5aa' for A, A' = 1,2,3. The first term in satisfies 
(□ + fi^)A^ = and d ■ A = 0. The A = term violates each, but the field equation is 
satisfied. The nonzero commutators are [a\{k), a\,{k/)] = e^5\\'S{k — k/). The resulting 
light-front Hamiltonian is 



dxTi 



x+=0 



dkJ2 



,^^±dalik)a,{k). 



(3.7) 



The nondynamical components of the fermion fields satisfy the constraints (i = 0, 1, 2) 



-l)^9_V'.-+eA_v^5^V^,-_ = (Z7V 



When multiplied by (— 1)*a/^ and summed over i, this becomes 



l)*mi7>i+. 
(3.8) 

(3.9) 
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which is the same constraint as for a free fermion field, in any gauge. 

Thus, the hght-front QED Hamiltonian can be constructed in any covariant gauge, to 
obtain 



+ j dydk± J 

ijlsaX 

xalxiyP'^^yP-L + k±)bl^{{l -y)p^, (1 - y)p± - k^)hi^{p) + H.c.j 



The functions h1^^ are readily computed from spinor matrix elements. Terms with positron 
contributions have been neglected, for simplicity. 



IV. APPLICATION OF THE LFCC METHOD 



We briefly discuss an application to the calculation of the anomalous magnetic moment 
of the electron. For simplicity, Fock space is truncated to exclude positrons. The valence 
state is |0^(P)) = Xli ^ib\±{P)\^)- The T operator is truncated to be 

T=Y.j ^^^^^ / -^^^Vp^t'i;i\y.kMx(.yp^^yp^+kMjs{{^-y)p^A^-^^^ 

ijlsaX 

(4.1) 

After exclusion of terms that annihilate the valence state, the effective Hamiltonian is 

,2 I ^2 



+ J dydk± J 



dp 



ijlsaX 

2 I U2 



+ 



mj + kj 1x1 + kj 

1 -y y 



hi;t{y,k^) + ^vr,t\yM) 



i' ^ j' ) 



xo-lxiyP'^^ yP^ + ki.)h]^{{l - y)p^, (1 - y)p± - k±)bia{p), 
with the self-energy Iji and the vertex loop correction V^'^i^ defined in ^ 



The projection P^V |0^(P)) 



'{P)) of the eigenvalue problem yields 



2 ± 



+ E ^^^4 = Ml 



(4.2) 



(4.3) 



Projection onto |e7), orthogonal to gives 



2 Mi + kl fil + k] 



i-y 



C^:,\y,k^) = Ht:i\yM)+ 



±sXi 



v^u\yM)-Y.Y^,Ca:'hyrk 



(4.4) 
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with 

C^u\y. ki.) ^ Y.i-iy^'z^rz^^tf^fiy, k^). (4.5) 

Notice that C^*'^, the analog of the two-body wave function, satisfies an equation where 
the bare masses have been replaced by the physical masses Ma. Thus, the LFCC method 
provides a natural way for physical masses to enter the calculation instead of having them 
imposed by a sector-dependent renormalization. Also, although the T operator has been 
truncated to a single photon emission, the projected eigenvalue problem retains the self- 
energy corrections and vertex corrections necessary for the Ward identity to be satisfied. 



V. SUMMARY 



Details of the application to QED can be found in [IJ]. The result for the anomalous 
moment agrees with the Schwinger correction at first order and will sum to all orders in a 
many of the higher-order corrections. 

The key advantage of the LFCC method is that it avoids Fock-space truncations that can 
induce uncanceled divergences and other inconsistencies, such as Fock-sector dependence 
and spectator dependence. The approach is systematically improvable, by the addition of 
terms to the truncated T operator. 

The complete calculation of the dressed-electron state and its anomalous moment, within 
the given truncation of T, is in progress. It requires a numerical solution of the right-hand 
and left-hand eigenvalue problems, as well as numerical quadrature for the matrix element 
that yields the anomalous moment. Work beyond this can include investigation of the 
dressed-photon state, extension of the dressed-electron state to include e+e~ pairs, and study 
of muonium, positronium, and symmetry breaking in scalar theories. Of particular interest, 
of course, is the adaptation of these methods to QCD; the LFCC method is general enough 
that it can be applied once a suitably regulated light-front Hamiltonian is constructed. 
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